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Abstract 

We apply the adiabatic approximation to investigate the low energy 
dynamics of vortices in the parity invariant double self-dual Higgs model 
with only mutual Chern-Simons interaction. When distances between 
solitons are large they are particles subject to the mutual interaction. 
The dual formulation of the model is derived to explain the sign of the 
statistical interaction. When vortices of different types pass one through 
another they behave like charged particles in magnetic field. They can 
form a bound state due to the mutual magnetic trapping. Vortices of the 
same type exhibit no statistical interaction. Their short range interactions 
are analysed. Possible quantum effects due to the finite width of vortices 
are discussed. 



Introduction 

Experiments with high temperature superconductors seem to show no indication 
of parity breaking Q . At the first sight this result seems to exclude the anyonic 
mechanism of superconductivity. But it is not the case as it was shown recently 
||. The presence of Chern-Simons interaction in a model does not lead 
inevitably to the breaking of the P and T invariance. It just a property of 
the simplest models with only one Chern-Simons field. When there is an even 
number of Chern-Simons fields and their coupling constants are appropriately 
choosen the parity invariance can be restored ^ ^| . One of the simplest models 
of this kind is the [[/(l)] 2 model of two Higgs fields each of them coupled to one 
of the two Chern-Simons fields Q. The model is constructed is such a way that 
particles which carry one kind of charge interact with the magnetic field of the 
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other kind. So the ordinary construction of anyons as charged particles which 
have at the same time attached a fictitious magnetic flux is split into to parts. 
Particles of one kind are carriers of the charge while those of the other kind 
carry the magnetic flux. Thus the ordinary fractional statistics (SJ is replaced 
by the so called mutual statistics 0. 

In this article we would like to investigate the Chern-Simons interactions 
of vortices in the relativistic self-dual [£7(1)]^ model Q. We apply Manton's 
idea of adiabatic approximation Q , with neccesary in the case of Chern-Simons 
vortices corrections Q], to the topological solitons of the model. The topological 
solitons configurations satisfy the lower Bogomolny bound on energy Q so the 
moduli space approximation is justified. Explicite calculations are made in 
the special case of the [[/(l)] 2 model but the results can be easily generalised 
to the case of [1/(1)]^. At large separations the vortices of different Higgs 
fields exhibit the expected mutual interactions but with a sign opposite to that 
expected from their fluxes and charges. It is very much like as for vortices in 
the standard U(l) Chern-Simons-Higgs model which exhibit ordinary fractional 
statistics properties ||, fjj. We derive dual formulation of the model to explain 
the sign of the statistical interaction. 

When the vortices pass one through another their interaction is a little more 
complicated. The pair of vortices of different types behave like charged particles 
passing the flux of magnetic field similarly as vortices in ordinary Chern-Simons- 
Higgs system j?], [|. Due to the fact that the spin of separated vortices is 
equal to zero while that of coincident is nonzero there is a kind of magnetic 
trapping - they form a composite state. If the corrections to the standard 
adiabatic approximation are quantitative in nature there is a periodic solution 
with vortices circling in the magnetic field of the trap. On the other hand the 
vortices of the same kind do not interact through the Chern-Simons field. If 
the corrections to the ordinary adiabatic approximation amount only to the 
rcnormalisation of parameters in effective Lagrangean, they would behave very 
much like vortices in the Abelian Higgs model. In particular the result of their 
head on collision would be the right angle scattering (h], [l2) . 

The article is organised as follows. In Section 1 we derive the general form 
of the effective Lagrangean. Section 2 is devoted to the long range interactions 
of vortices and their mutual statistics. In the next paragraph (3) we analyse 
what happens when various types of vortices pass one through another. Section 
4 is a presentation of the dual formulation of the model. In the last section we 
summarise and discuss the results. 

1 General form of the effective Lagrangean 

The Lagrangean of the relativistic model presented in Q when we restrict to 
such a choice of parameters that only mutual interactions are preserved can be 
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written in the form: 



L = KE^ x v^d v vf ] + D^*D^ + D^D^x - V{4>, x) , (1) 
where the covariant derivatives are defined by 

D M <j> = d^-iqiv^cj) and D^x = d^X - mv^x > ( 2 ) 
while the Higgs potential of this self-dual model is equal to 

Vfax) = #11 I 2 (I X I 2 ~4f+ | X I 2 (I 4> I 2 ~c\f] . (3) 

First we will work out the general form of the effective Lagrangean by a direct 
application of the methods of adiabatic approximation known from the papers 
on slow motion of vortices in self-dual Maxwell-Higgs system 12 and then 
we will discuss corrections to this oversimplified version of the Lagrangean . 

The Lagrangean (|l|), when we eliminate auxiliary components of the 
gauge fields with the help of the Gauss law 



(i) q kB 2 ( 2 ) „ kB-i 

2qi \0\ 2 2g 2 x 1 



can be rewritten in the form 

L = -ns^v^dtv^ + —B (2) d t uJi + —B (1) d t u>2 
J qi 92 

+ (dt\0\) 2 + (d t \ X \) 2 

-i^-i^-\^\ 2 -\Dx\ 2 ~V^x) • (5) 

where we have introduced magnetic fields corresponding to the two different 
gauge fields: Bj = e^diVj. wr-s are the phases of the Higgs fields. We have 
separated terms containing time derivatives from those with spatial derivatives 
of the fields. In this form of the Lagrangean we can make replacements due to 
the identity 

- | | 2 - - | ( J D 1 +i<T 1 £» 2 )0 | 2 -a iqi B w | 4> | 2 +ia 1 s ij d i (ct>*D j <j>- c.c. ) (6) 

and an analogous one for the field x- The constants a\ 1 <j 2 which take values 
from the set { + 1,-1} have the same sign as the topological indices of field 
configurations cj) and x respectively. 

In the Manton's approximation we assume the fields to have the same form 
as in some static solution and thus to fulfil the static field equations at every 
instant of time. In our case the fields fulfil the self-dual equations: 

(Di + i<TiD 2 )<j) = and (D 1 + ia 2 D 2 ) X = (7) 
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so the first term on the R.H.S. of eq.(|^) vanishes. The assumed static solution 
depends in general on some finite set of parameters which will play the role of 
collective coordinates. To derive an effective Lagrangean we have to integrate 
out the spatial dependence of the fields in eq.(0). Assuming the Coulomb gauge 
condition for the static solutions, which implies that the gauge potentials can 
be written as 

v\ r >=e ii d i vn (8) 

with some regular functions we can remove first term in the first line of 

(||) and third term on RHS of Eq.(||). Further simplification can be achieved 
with a help of the remaining static field equations 

B (1) = -2a 1 ^| X | 2 (U| 2 -c?) , (9) 

B m = -2a 2 ^|l | 4>\ 2 (I X I 2 -cl) • (10) 

Neglecting terms which are constant in a given topological sector we obtain the 
following form of the effective Lagrangean 

L ef f = I d 2 x p t | |) 2 + (d t | X I) 2 + -BadtU! + -B 1 d t u 2 ] . (11) 
13 J qi Q2 

Y.Kim and K.Lce M tought us that it is important to take into account cor- 
rections to such direct application of adiabatic approximation which is oversim- 
plified in the case of Chern-Simons vortices. Namely even in the slow motion 
approximation the fields can not be taken just as static solutions with time 
dependent parameters. We also have to take into account corrections to those 
fields which are linear in velocities. 



X 



X + S X , v<P^vjp+8vP (12) 



because they give additional terms to the effective Lagrangean. It is a kind of 
complicated Lorenz transformation. 

We can show that corrections affect only this part of effective Lagrangean 
which is quadratic in velocities. Let us take closer look at the Lagrangean 
in the form of Eq.(||). Second line of the above formula is manifestly of at 
least second order in velocities. Third line is minus static energy density, so 
the static solutions are its stationary points. Thus third line also gives only 
quadratic terms. The only contribution to the linear part can come from the 



first line. Any replacements of the form (12) will produce only extra second or 
higher order terms. So we have to take the first line as it stands - only with 
static fields with time-dependent parameters. As it has been shown the first 
term of the first line gives no contribution, so the linear part of the Lagrangean 
remains as 

L« = K f fa^O-fa + ^ W2 ] . (13) 
JJ J qi 92 
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Another representation of L^jf, useful for further discussion can be obtained 
with use of the identity 

ui = 4; Arg(x - Rp,) = 01 e^R % m dj In [ x - Rp, | (14) 

pi pi 

and certain integration by parts 

4 1 / ) / = ^E^f ) (^)+ (15) 
1 pi 

R Pl (t) are positions of vortices of type "I" ,1=1,2. From this representation one 
can obtain general form of the orbital part of angular momentum 

~ Z^Z^^^Tl an) 
I pi ° n Pi 

2-KKp-Y. R v^t\K)+ (1-2)] . (16) 
91 pi 

The form of L^. h can be extracted from (||) after replacement (|l^) and making 
use of static equations fulfiled by "static" fields. It gives correction to the 
effective Lagrangean ( p"T| ) 

5L e f f [fields, S. fields] (17) 

which is a functional of "static" fields, their linear corrections and their time 
derivatives. It is a long expression which we will not write down. To make 
use of this expression we have to express "5-fields" in terms of "fields" and/or 
positions of vortices. We have to know corrections for a given trajectory in the 
parameter space. We can make replacement ([l2]) in the field equations following 
from (|) 

D M £>^=-^ , Ke^ x d v v ( ? =%(0D^*-c.c.) (18) 

and analogous one for the field X- It is convenient to use geodesic parametri- 
sation on moduli space. By geodesic parametrisation we mean such a set of 
parameters that during the time evolution of the field configuration their ac- 
celerations are vanishing. It is always in principle possible to construct such 
geodesic coordinates at least for a finite period of time. One advantage of such 



a set of coordinates is that we can neglect in Eqs.(18) terms which contain accel- 
erations or higher order time derivatives of parameters. The other one also very 
important is that only when velocities in a given parametrisation are constant 
there is a direct correspondence between the fact that kinetic energy is small and 
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smallness of velocities. For example in the head-on collision of Nielsen-Olescn 
vortices the time derivatives of cartesian coordinates of vortices become singular 
during collision |kJ 12 1, so in this case terms quadratic in such velocities are 



much larger then linear terms and obviously can not be neglected. Nethertheles 
the kinetic energy is still small - it is still a slow motion. In fact it is possible to 
find geodesic coordinates good for small separations of these vortices. So only 
for geodesic parametrisation it is reliable to preserve in Eqs.([l8|) only terms 
which are linear in velocities and neglect accelerations and higher order time 
derivatives. 

Thus "simplified" equations read 



= d 2 k 5h + - (d^ - qiv^f]5h + 2q 1 f 1 (d k cj 1 - qiv^)6v^ 

d 2 V d 2 V 
dfl h dhdh' 



-2?i/i^ ; (w ( i) - giSv y Q ') - jr^Sh - 8 f a 



= fM 1] +d k Svf 1) ) + 2q 1 v ( ^f 1 

= KEijdMp - qifl^i - qiSv^) + 2q 2 1 v£ ) f 1 5f 1 

= KEijidjSvjP - vf ] ) + qlfitv^ -2q 1 f 1 (d k Lu 1 -q 1 vi 1) )6f 1 (19) 

and an analogous set of equations for the field x- We have introduced moduli 
and phases of the Higgs fields: = -^e 4 " 1 , X = "^ e ^ 2 - ^ n analogous set of 
equations was derived in [Q for self-dual Chern-Simons-Higgs model. Neglecting 
of second order time derivatives for any kind of coordinates can lead to serious 
problems. For example in the papers on effective string models for Nielsen- 
Olesen vortices the second order derivatives on world-sheet parameters were 
neglected. It was shown in Jl6| that it was the reason why the string with 
rigidity was obtained, which is known to posses classical tachionic solutions 

With a help of Eqs.(|l|) SL e ff can be simplified to the following form 
i J 

-fiSftidiU! - qivP) qi SvP (20) 

Once we have solved Eqs . dl9|) we can substitute their solutions to the above 
functional and integrate out their spatial dependence. We will be left with a 
mechanical Lagrangean which should be reparametrisation invariant. From that 
point on we will be able to use any parametrisation we like. But to evaluate 
SL e ff we have first to guess geodesic coordinates and then to solve Eqs.([l9|). 
We have to try with different parametrisations and then to check whether for 
a given parametrisation there exist solution to Eqs([l9|). If yes than as a matter 
of fact we have found already the trajectory. The effective Lagrangean can be 
evaluated for consistency check and because it is useful if we want to perform 



G 



effective quantisation of the theory or to investigate its effective thermodynamics 
p3). The effective Lagrangean under restrictions due to the above comments 
reads 




To start the above procedure we have to make some "educated" guesswork. 
To provide an appropriate basis for it in the next two sections we investigate 
ordinary adiabatic approximation in detail. Because of mathematical difficulties 
(existence proof and/or explicite solution of Eqs.([j"9|)) we postpone calculation 
of corrections to future publication. 



2 Long range interactions of solitons 

For sufficiently separated vortices we can approximate gauge invariant fields: 
(dkUJi — qivl) by contributions due to particular vortices. At the core of any 
choosen vortex such fields due to the other vortices are very small - they vanish 
exponentially with distances. Since u>i — ai ^2 PI Arg(x — R Pl ) is an exact 
solution of static equations for any configuration of vortices, corrections to the 
simple superposition of gauge potentials are negligible. Thus values of the gauge 
fields in Eq.([l5|) can be obtained from a formula 

vf\R Vl ) = -^Y j £ij{ ^-^ )i , (1-2) (22) 
so the linear part of effective Lagrangean is 

Pl = lf»2 = l 

The Lagrangean contains only terms of mutual statistical interaction between 
vortices of different types. 

To obtain the kinetic term one has to evaluate the first term of the effective 
lagrangean (|2l]) and SL e f j . We approximate the moduli of the Higgs fields by 
a normalised product of the fields of isolated unit vortices: 

m 

\4>{x)\=ciT{G{\Z-R$\) , (24) 

Pl = l 
712 

\ X (Z)\=c 2 ]jG(\x-R$\) , (25) 

P2 = l 
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where G is a profile of a unit vortex the same for the two types, which satisfies 
the equation 

V 2 In G 2 = M£l£M (G 2 - 1) , (26) 

K 

with the boundary conditions: G(0) = and G(oo) = 1. The quadratic part of 
the effective Lagrangean reads 

L % - 2^ E i£ + tg] + «W/ , (27) 

Pl = l P2 = l 

where the coefficient M equals to 

M = 2n[rdr(^p-f . (28) 

About SL e ff we know only that it is quadratic in velocities. Let us restrict to 
the case of two vortices of the same type and choose the center of mass frame 

i.; fi u. Air .iriv , 9ij = 9ji , (29) 

where Rk-s are coordinates of the choosen vortex. Rotational invariance restricts 
this form to 

L^ f = 9l (R)R 2 +g 2 (R)R 2 Q 2 , (30) 

where R 1 + iR 2 = Re lQ . At very large R we expect the influence of one vortex 
on another to be very small: gi(R) , 92{R) — > const as R — > oo. This reasoning 
can be repeated for any pair of vortices. Finally we obtain 

Leff = \m\c\ £ Vl +c 2 2 J2 V 2 2 ] + Yl E J t A ^ - 3g>) ■ 

pi = l P2 = l Pl=lP2=l 

(31) 

where M = M + JM is the effective mass with included corrections from SL e f f . 
Corrections can be calculated with a help of the formula (^), since for fairly 
separated vortices " 6- fields" are given by ordinary Lorenz formulas linearised in 
velocities. The corrected coefficient M appears to be equal 2ir, what is consistent 
with the original field theoretical model since the energy of a static unit vortex 
of type "1" is equal to 27rc 2 . 

We can see that when the widths of the vortices can be neglected as compared 
with their separations the system behaves like a set of free particles with mutual 
statistical interactions, at least in the slow motion approximation. What hap- 
pens if the vortices come into very close encounters of one another is a subject 
of the next section. 

Let us remark here on the possibility of ordinary fractional statistics in the 
system if the short range interactions between vortices of different types fauvored 
them to form mixed anyonic cj> — \ bound states. 
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3 Short range interactions 



In this section we would like to investigate interactions of the two types of 
vortices when their cores overlap. To apply the Manton's prescription we have 
to know at least an approximate static solution. Let us take as a zero order 
approximation the configuration of vortices sitting on top of each other and 
then let us find a small perturbation of the field dependent on a definite set of 
parameters. 

Let the topological indices of the fields <f> and \ be <7\n\ and <72"-2 recpec- 
tively, where n-s are positive integers and ct-s take values +1 or —1. The solution 
corresponding to vortices sitting on top of each other takes the form: 

<j> = Cl F(r)e l,7inie , X = c 2 H{r)e ta2n2e , (32) 

"(l) = eeCTiV^) , V( 2 ) = e e a 2 V® 2) . (33) 

Upon substitution of the above Ansatz to the static field equations (0,^,0) one 
obtains 

F' — —F + qiFVg 1 ^ =0 , (34) 
H'-^H + q 2 HV g (2) = , (35) 
-ld(rV g (1) ) 



q 2l H 2 (F 2 -l) , (36) 



r dr 

±2£2&_ B7na ._ 1)l7 _2£«££i. (37, 

r or k 
This zero order solution together with a small perturbation reads 

(j) = Cl F(r)[l + f(r,8)]e Unnie+iai W , (38) 

X = c 2 H{r)[l + h{r,9)}e la2n2e+la2{r ' e) , (39) 

v {1) = eeaiVfa + a (z) , I = 1,2 (40) 

Upon substitution of the above form of the solution to the self-dual equations 
(0) and subsequent linearisation one obtains the following first order equations 

qja^ = U e a^ - ajd r f^ , (41) 

qi ai I )=d r a^+a I ^- . (42) 
r 

These equations enable us to find perturbations of the gauge fields once = 
(/, h) and ot(iys are already known. 
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After substitution of eqs.(|4C) to the Coulomb gauge condition, diV^ 1 ' = 0, 
linearisation and use of eqs.(|4l|^42|) an equation determining the phases of the 
Higgs fields appears: 

V 2 a (/) = . (43) 
Similarly equations ([34 35 ) linearised in the perturbations yield 



V 2 / =jH 2 [F 2 f+ (F 2 - l)h] , (44) 

W 2 h = jF 2 [H 2 h+{H 2 -l)f] , (45) 

where 7 = 2q\q2 r y. From the two above equations one can obtain general form of 
the pertubation of the moduli of the Higgs fields. Then one have to choose such 
a solution of eq.(43) as to avoid singularities of the gauge fields, see eqs.( [4l]Jl2"| ). 



We will solve this problem explicitely in two special cases. 

3.1 Interaction of unit 0- vortex with unit \~ vortex. 

In this case we have unit topological indices of both the </>-field configuration 
and the x-field one 

ni = 1 , n 2 = 1 . (46) 
Eqs.gl^) can be simplified by a substitution 

F(r) = H(r) = G(r) , S- = -J- = V e (r) , (47) 

to the form 

£?'- — + qiqvGVg = , (48) 
r 

Z 1^) =7G2(G 2_ 1) (4g) 

r or 

From the index theorem Q we know that in the special case of ii\ = 712 = 1 
there are only two splitting modes for each of the two types of the fields. We 
can Fourier-transform / and h in 0: 

f(r, 0) = f{r)[\ x cos{ax0) + X 2 stn(a 1 9)] , (50) 

h(r, 9) = h(r)[/iicos(a 2 0) + ^2sin(a 2 9)} (51) 

The ^-independent terms are neglected because solitons have a definite size, 
while the terms higher than the first would spoil regularity of the total Higgs 
fields. To define the meaning of coefficiens A, fi we normalise the radial functions 
in such a way that 

f(r) ~ — , h(r) ~ — ,asr^0 . (52) 
r r 
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Eventual higher order terms would have a stronger singularity which could not 
be matched by G(r) ~ r, see eqs.( |38| , |39| ). Thus eqs.( |5l] , p2| ) present the general 
form of the perturbation compatible with the regularity of Higgs fields. 

To avoid singularities in the gauge fields, eqs.(|4l],^), we take the perturba- 
tions of phases 

ct\ = — \\2Cos(ai9) — Aism(<72#)l , (53) 
r 

«2 = — [fi2cos(a 2 0) - /iism(a- 2 0)] , (54) 
r 



which satisfy eqs.((43[). Now we can see that the perturbed Higgs fields in the 
limit of small r are proportional to 

<f> ~ (z ai - A) , x - (z a2 - M) , (55) 

where we have introduced z ai — x + iajy and A = Ai + 1X2 , (i = /ii + 1^2- 
The effect of the perturbation is a shift of the zeros of the Higgs fields to A CTl 
and to /v 2 1 U P to linear terms. To work in the center of mass frame we have to 
choose /v 2 = — A (Tl together with c\ — c-i — c. Last condition can be suspended 
if we make certain relative rescaling of the parameters A and /i. Without loose 
of generality in evaluation of f(r) and h(r) we can choose X2 = ^2 = 0. Now 
upon substitution of eqs. (|50| , |5T| ) to eqs.(^J^) one obtains 

A 1 /(r)=7G 2 [G 2 /W-(G 2 -l)/ i M] , (56) 

A 1 h(r) = jG 2 [G 2 h(r) (G 2 - l)/(r)] , A fc = (-J + ±A - %) . (57) 

ar z r ar r z 

We can make the replacements f(r) = w(r) + u(r) and h(r) — w{r) — u{r), 
where the newly introduced functions fulfil the equations 

Aiw(r) = jG 2 w(r) , (58) 

A lU (r) = jG 2 (2G 2 - l)u(r) . (59) 

The only solution compatible with normalisation ( p^ ) is u(r) = and f(r) — 
h(r) = w(r). w(r) can be approximated for large r by the modified Bessel 
function K\{p) ~ ^J^ e ^ p >P = \fl r - Going from infinity to zero this approx- 
imate solution is replaced by a linear combination 5+1?" + &—\-. In the case of 
7^ 0, which we think to be quite general, it is possible to rescale the whole 
function w(r) in such a way that we obtain an assymptotics of eq.(|52|). The 
perturbations of the Higgs fields are square-integrable. 
Now the effective Lagrangean @ reads 

Leff = \M eff {R 2 + R 2 Q 2 ) + 8L eff + B eff R 2 <3 . (60) 
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R and are the polar coordinates of the zero of the 0- field: A = Re , while 
the introduced coefficients are an effective reduced mass of the two vortices 

/>OC 

M e// =47rc 2 / (rdr)G 2 (r)w 2 (r) , (61) 
Jo 

and an effective "uniform external magnetic field" 

B« t =*^r*&(r)[-w(r)] . (62) 
Kaia 2 Jo 

Let us first discuss linear part of the Lagrangean (|60|), which becomes exact for 
very small R. It is a term which describes, as it stands, coupling of a charged 
particle to a uniform external magnetic field perpendicular to the plane. The 
total angular momentum in the effective description is (for = 0, R — ► 0): 

J(R) d g 2« + dL^ = _2nna 1 a l + ^ + 0{r3) 

where we have shifted the scale so that for R — > 0, J tends to the value of spin 
characteristic for coincident static unit <fi and x vortices. With this choice of the 
scale we obtain from (|3l|) that for large R: J — > 0. This result is consistent with 
what we know from field-theoretical considerations. Separate <fi or x mutually 
interacting vortices carry no spin. Spin is nonzero only when their cores overlap. 
Eq. ( |63| ) gives leading terms in expansion of J{R) around R = 0. 
It is interesting that Eq.(|63|) can be inverted in a remarkable way 

1 ,J(R) - J(0), ,. 1 dJ(R) 
B eff = lm — (-^ — ) = lim ^ . 64) 

en r^o R K R ' n.->o R dR y ' 



Natural thing is to ask whether such a formula can be generalised to arbitrary 
value of R. Let us look at in the form of Eq.([l5|) 

L « = ^^^f)(4 i)+ (1++2) . ( 6 5) 
1 pi 



What we see is an interaction term which couples point particle currents to 
fields defin 
magnetic field 



fields defined on the moduli space. Due to this coupling vortex at R pi feels 



B« (R P1 ) = 2 -^l £ij d i vf\R Pl ) = (66) 

Our pair of vortices feels double this field: -B e // = ^e// + ^eff' From the 
formula ( |l6] ) we obtain angular momentum, which in our case reads 

J orb (R) = 2%K—R i E i jvf ) {R)-2TrK—R i ei j v\ 1) (-R) , (67) 
<7i <?2 
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Due to the rotational symmetry, in polar coordinates 



Jorb{R) = 2-KK—Rvf\R) + 2ttk—Rv£ ) {R) . (68) 

Now we can see that 

1 dJ 1 dJ orb oi ( vf ] , dvf\ 

= = 2ttk — (— — + ( 1 •*-> 2 ) = 

RdR R dR qi R dR ' V ' 

= fl« (R) + B% (R) = B eff (R) (69) 
Thus eq. (|64|) can indeed be generalised for any value of R: 

1 dJ(R) 

Beff{R) ~R^R- ' (?0) 

where B e ff{R) is a magnetic field felt by a reduced particle in our problem, 
while J(R) is a total field-theoretical spin, which can be obtained numerically. 
It is important that this formula is based on global properties of the field config- 
urations and not on local distortions of the fields, so B e ff(R) can be calculated 
numerically with great accuracy. 

Just from the knowledge of spin dependence on R we can obtain qualita- 
tive understanding of interactions between <f> and x vortices. Vortices boosted 
against each other for a head-on collision will avoid direct collision. Their initial 
total angular momentum is zero. For very small separations spin itself would 
have to be close to (— 2 ^^"^ 2 ) so the vortices must aquire angular momentum 
of the opposite value. For 2^2. positive they will be turned by B e ff(R) to the 
right while for the negative value of the coefficient they will turn left. Since 
B e ff(R) is short ranged, for sufficiently small initial velocity it becomes impos- 
sible to reach the center of mass for any value of impact parameter. Analogous 
argument shows that if a vortex initially sits on top of \ vortex there is a 
velocity small enough below which they can not escape to infinity. So it is a 
kind of magnetic trap for vortices. They trap each other and form a compos- 



ite. From Eq.(31) is clear that such composites as a whole behave like anyons. 
Their internal reduced dynamics is that of a particle in an external magnetic 
field. Upon quantisation it is probable to obtain resonances classified by the 
Landau levels. The trapped vortices do not form a real bound state - their 
energy is the same as that of isolated vortices. 

Now let us take a look at the quadratic part of the effective Lagrangean. If 
we assume that SL e f / only renormalises M e f / than L e f / is a Lagrangean for a 
planar motion of charged particle in external uniform magnetic field perpendic- 
ular to the plane. Solutions to the equations of motion are circular trajectories 

R(t)e i@ W = A + Be- mt , SI = , (71) 

M eff 
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where A, B are complex constants. This solution is valid only for small R. 
The trapped vortices rotate around circles (not necessarily around the center of 
mass if A ^ 0). We can also obtain qualitative understanding of scattering. Let 
us assume that this approximation is valid up to say i?o ; for larger R let the 
vortices move along straight lines. In the head-on collision vortices cross the 
circle R = Rq, move along an arc of a circle and escape to infinity. For larger 
velocities, if the adiabatic approximation still works, scattering pattern evolves 
to a forward scattering. We can check whether M e /y is indeed a constant for 
small R by performing numerical simulation of the above described head-on 
collision. For example if for higher velocities right-angle scattering is obtained, 
M e f f must behave like R 2 for small R. From analogous simulations for Nielsen- 
Olesen vortices we can expect that the adiabatic approximation can still work 
well even up to | of the light velocity. 

The eventual dependence of M e ff on R or its constant renormalised value 
can be obtained from a measurement of the frequency of the purely orbital 
motion of the trapped <fi ~ X pah: 

^'<*» = W' ,72) 



where it is assumed that B e ff(R) is already known from ( |70| ) and fl(R) has 
to be measured in numerical simulation. Alternative analytical approach is to 
take O as a geodesic coordinate and try to find such a value of f2 = for which 



Eqs.(19) posses a uniquae solution. 



3.2 Interactions of unit vortices of a given type. 

It this subsection we would like to investigate short range interactions of say </>- 
vortices when there are no x~ vortices or their influence can be neglected because 
they are very distant. We put \ — c i and n\ — n , a\ — +1 in this paragraph. 
The configuration of n vortices splitting from their coincident position is of the 
form 

rtr,9)=cF(r){l + f(r,8)}e me+la ^ , =V= e g V e (r) + a(r, 6) , (73) 

where the functions F(r) and V (r) satisfy equations 

ft 

F' F + qi FV = , (74) 

r 

-ld(rV e ) 2, P 2 ,s 

~ 1ilc 2 {F - 1) , (75) 



r dr 

Following similar steps as in subsection 3.1 we obtain equations fulfilled by the 
perturbations: 

q 1 a i = d i a + Sijdjf , (76) 
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V 2 a = , (77) 

V 2 / - iF A f (78) 

and solve them by an Ansatz 

f(r, 9) = g(r)[X lC os(n9) + X 2 sin{n9)} , (79) 

a(r,9) = ^-[X 2 cos(n9) - X 1 sm(n9)} , (80) 

where g(r) is a solution of equation 

A n g(r)=iF i (r)g(r) . (81) 

g(r) is normalised so that g(r) ~ ^ for small r. We have taken the n-th term 
in the Fourier transform because it corresponds to a uniform spliting of the n 
vortices from a coincident position 

4>~ (z n — X) , A = Ai + iX 2 . (82) 

The positions of the n vortices are n-th order roots of A, which we denote 
by Re l ® +l ~^~ k ,k — 0, (n — 1). In the case of only 0-vortices the effective 
Lagrangean reduces to 

L eff = J d 2 x \d t ^\ 2 +SL eff =ncj J rdr F 2 g 2 {X 2 +X 2 ) + SL eff . (83) 

When we take into account that A = R n e me we can rewrite the Lagrangean as 

Leff = ~M e// (AA*) + SL eff = ^M eff R 2 ^(R 2 + R 2 Q 2 ) + 5L eff . (84) 

If we neglected SL e f / in (§)Jt would follow that in the head-on collision A turns 
to —A and it is clear from (|8^) that it means scattering by an angle of — . The 
configuration of n vortices shrinks to a coicident position and than reappears 
but rotated by an angle of ^ with respect to the initial one. 

Let us analyse corrections due to SL e ff. First we make a hypothesis that 
A is a geodesic coordinate or A = for values of A. The "fields" are regular as 
functions of A in A = 0. By regularity we mean finiteness and single- valuedness. 
The " (5-fields" are defined as linear in A and we take them as a series in the 
powers of A (components Ai,A2). If we substitute such "fields" and "<5-fields" 
to the field equations ( |l8|) we will find that because "fields" are regular the 
"(5-fields" can also be taken self-consistently as regular in A = 0. In the limit 
A — > we will obtain equations which are linear in Ai,A2 and in the values of 
"(5-fields" at A = 0. If there is a solution to these equations we can use it to 
calculate 5L e ff and because "(5-fields" are regular at A = correction to the 
effective Lagrangean amounts only to renormalisation of M e ff. If there is no 
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solution it will mean that the initial assumption A = was wrong and we have 
to look for other candidates for geodesic coordinates. 

Following argument can restrict the set of acceptable candidates. Let us 
take the head-on collision of two 0-vortices. Initially they are coming to the 
center of mass along x-axis. Such a configuration is invariant under succesive 
charge conjugation and reflection with respect to the y-axis. Since our theory 
is CP-invariant the time evolution has to preserve this symmetry of the initial 
configuration, so if the zeros of the Higgs field pass through the center of mass 
there is possible only forward scattering or right-angle scattering. For the right- 
angle scattering A,-s are good geodesic coordinates but the forward scattering is 
well described by the cartesian coordinates of vortices (zeros of the Higgs field) . 
If A = leads to contradiction that means that we will have to try with R, O. 

Since these calculations are a fairly nontrivial problem we will only conclude 
that only forward or - scattering are possible in the head-on collision of vortices 
of the same type. Knowledge of the exact form of L e tf could be useful to 
effective quantisation of the model and to investigations of its thermodynamics 




4 Dual formulation 

In this section we derive dual formulation for mutually interacting vortices, 
following similar steps as in for ordinary Chcrn-Simons-Higgs system. There 
are two reasons for it. We would like to show that the dual transformation 
can be in a natural way generalised to the systems with more complicated 
Chern-Simons terms. Second and more important - it explains why the sign of 
statistical interaction is inverse to what could be expected from calculation of 
naive Aharonov-Bohm phase. We will use Lagrangean ([[]) with extra couplings 
to external currents and external field 

L = KE^^duV^ -V{f X ,h) + 

D^^^i-^-^ + Wi)] > ( 85 ) 
i=i 

where A e * is an external field and are external currents. We have rewritten 
Higgs fields as 

0=-L/ lC <Wl , x = -L/ 2 e^ . (86) 
The partition function is 



Z = 



J JJ[4fj] [dui\ W;\ \{ exp iJd 3 xL . (87) 



16 



Phases of the Higgs fields can be split into multivalued and regular parts 

u)i(x) = u>i(x) +r)i(x) , (88) 

where Q-s are given by 

u I {t,r) = Y J ° Pl Arg{r-R Pl (t)) ■ (89) 
pi 

R PI are positions of unit vortices (a Pl — 1) and unit antivortices (o Pl = — 1) of 
the type " I" . It is covenient to construct out of u> vortex currents 

i r dR^ (t) 

K$(x) = —e» vX d v d x m = $>« / dT P d ' T < x ~ R p'^ ( 9 °) 

pi 

in a covariant fasion. By definition K satisfies the conservation law: d^K^ = . 
Integration over ui can be replaced by integration over vortex worldlines 

[duj} = ldu>}{dr,} = {dR£ I }{d V } • (91) 

The Jacobian F| 7 FT // (x) can be removed by introducing pair of auxiliary fields 
Of: 



nil expifcPx^fHd^ - qi vl - ei Aff = 

/ X 

J XlMd] exp i J d^E/h^^C? + CWnQi + d^j - qjv^ - ejA^)} (92) 

Integration over 77-s will introduce Y\i 5{d^C^) to the path integral measure. 
These (^-functions can be removed by introducing a pair of dual gauge fields 

\[[dCf]5{d^) = J l[[dC?][dH?]6(C¥ ~ ^L-^duHl) (93) 

and integrating over auxiliary fields Now the vortex currents can be in- 

troduced by the identity 

[ d 3 x J— S ^d nQi&Hl = - [d 3 x Kf H 1 , (94) 
J 2irqi qi J 

where integration by parts has been done and we have made use of the definition 
of Kj , see ( po| ) . The present intermediate form of the Lagrangean reads 

1 

~V{h,f 2 ) (95) 
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and the integration measure 



nMf/«jwp<] • (96) 

I 

We would like to remove the gauge fields Vj from the Lagrangean. Let us take 
a look at their classical field equations following from (p5[): 

Ke^duvP + ^ vX d v H { ^ + J* = , (97) 

Z7T v ' 

Ke^ x d v v { ^ + —e^dvH^ + Jfa = , (98) 
Motivated by these equations we write 

»P = + > «£° = " a^4 l) + G f > CM) 

where GjP are extra fields due to the presence of external currents and contain- 
ing quantum fluctuations around classical solution. Thus finally the partition 
function reads 

Z = J YlWi^d^ldH^ldG^} exp i J d 3 x L D , (100) 
where the dual Lagrangean is 

+ S^(»»W , ]-f(/i,H] 

I 

+ne^ x G^ ) d u G^ ) 

D-^^W + ^ - i< J « - i< )j (^ 101) 

When there is no external current the fields G 1 ^ decouple and can be integrated 



out giving contribution to the normalisation factor. If in addition = wc 
are left only with the first two lines of the above dual Lagrangean. 

Vortex current couple to the dual fields . Since the Chern-Simons term 
for the dual fields has an opposite sign to that for the original gauge fields Vj it 
explains why the statistical interaction term in the effective Lagrangean (Q) has 
an opposite sign to that expected from the values of vortex fluxes and charges. 
The dual Aharonov-Bohm interaction between vortices is mediated by the dual 
fields and it gives rise to the correct value of the statistical interaction. 
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5 Conclusions 



We have made an analysis of interactions of self-dual Chern-Simons vortices in 
the limit of very large and very small separations. We have shown the existence 
of mutual statistical interaction between vortices of different types in [£/(l)] 2 
model but the results can be easily generalised to the general [(7(1)]^ theory. 
The sign of the statistical interaction is inverse to expectations based on ordinary 
Aharonov-Bohm effect. Thats why we have derived dual formulation of the 
system in which it is clear that vortices interact via dual gauge field with the 
sign of the mutual Chern-Simons term inverse to that in the original formulation. 

We have not attempted calculating corrections to the standard adiabatic 
approximation but a possible method how it could be done was discussed. If the 
corrections are only quantitative in nature the qualitative picture of short range 
interactions obtained in ordinary adiabatic approximation remains unchanged. 
In the head on collision of vortices of the same type we should observe right-angle 
scattering. For vortices of different types at large separations dual Aharonov- 
Bohm effect is observed but when their cores overlap they behave like charged 
particles crossing magnetic flux. 

The analysis of the short range interactions of vortices in Chcrn-Simons- 
Higgs systems presented in both this article and in |7|, |) shows the possibility 
of periodic solutions very much like bound states of vortices. The semiclassical 
quantisation of these solutions |ll[] can give rise to some discrete spectra of 
energy. The spectra can be expected as an additional quantum effect to the 
statistical interaction due to the finite width of vortices. The special effects of 
the short range interactions could be expected to vanish in the limit of vanishing 
thickness of vortices but we know from the studies of the string limit for vortices 
in Abelian Higgs model jl4j that even when classical vortices become very thin 
the quantum fluctuations cause that they preserve nonzero effective thickness. 
So it is possible that a classical vortex of finite width is a better zero order 
approximation to the full quantum theory fl5|| . 

The other topic worth of detailed study is the possibility of existence of 
ordinary fractional statistics in the apparently only mutually interacting system. 
The pair of ip and x vortices can form a composite thanks to the magnetic 
trapping. If the potential of the model V(4>, \) were slightly deformed in such 
a way that it would prefer energetically overlapping of the 0-vortices and %- 
vortices but it would discourage vortices of the same type to overlap, than we 
would expect vortices of different types to form true stable bound states. If we 
wanted this bound states to be composed of exactly one vortex of each type we 
would have to make a repulsion of the species of the same kind to be stronger 
than attraction of vortices of different types. Such a multivortex system in a 
sufficiently low temperature would be a gas of such anyonic bound states. In 
a higher temperature the average kinetic energy of the anyons could be large 
enough to split them into particular mutually interacting vortices. Thus we can 
construct a system with two phases: an anyonic one and a phase with mutual 
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statistics. The composed anyons might have an interesting internal structure. 
If the corrections to the Higgs potential do not change to much the interactions 
patterns at small separations the two vortices will feel both the charge-flux 
interaction and an oscilatory type interaction due to the atractive properties 
of the Higgs potential. The phase-diagram of the system could be even more 
complicated if the Higgs potential itself depended on temperature. We think all 
these topics to be worth of further investigation. 
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